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We are using domain wall fermions to study K — > tttv matrix elements by measuring K —* n and K — > matrix 
elements on the lattice and employing chiral perturbation theory to relate these to the desired physical result. The 
residual chiral symmetry breaking of domain wall fermions with a finite extent in the fifth dimension impacts these 
measurements. Using the Ward-Takahashi identities, we investigate residual chiral symmetry breaking effects for 
divergent quantities and study pathologies of the quenched approximation for small quark mass. We then discuss 
the AS = 1 operator O2, where chiral symmetry is vital for the subtraction of unphysical effects. 



1. INTRODUCTION 

Lattice QCD allows a first principles deter- 
mination of hadronic properties, including the 
matrix elements of operators between hadronic 
states. Among the more interesting matrix ele- 
ments are those of the AS = 1, four-fermion oper- 
ators which are responsible for the decay of a kaon 
into two pions. These matrix elements are needed 
to relate the values of the Cabbibo-Kobayashi- 
Maskawa matrix in the standard model to physi- 
cal measurements. Direct measurements of K — > 
7T7T on the lattice are not possible jlj, due to 
the two pions in the final state, although a re- 
cent suggestion involving tuning the lattice vol- 
ume Q looks promising. Alternatively, lowest 
order chiral perturbation theory can be used to 
relate K — > 7T7t to K — » 7r and K — > matrix ele- 
ments ||, which can be measured on the lattice. 
This approach has been extensively studied and, 
while quite straight-forward in principle, requires 
careful control over a wide range of phenomena. 

Of central importance to a calculation relying 
on chiral perturbation theory is precise control 
over chiral symmetry on the lattice. A major ad- 
vance in this long-standing problem is the idea 
of four-dimensional chiral fermions arising from 
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mass defects in a higher-dimensional fermion for- 
mulation With improved chiral symmetry at 
finite lattice spacing, the lattice Ward-Takahashi 
identities, without complicated renormalization, 
are much closer to the continuum ones. Here we 
seek to understand how the residual chiral sym- 
metry breaking effects for domain wall fermions 
enter the removal of unphysical, power divergent 
contributions, which is needed to get K — » 7T7t 
from K — > 7r and K — > matrix elements. 

2. DOMAIN WALL FERMIONS 

The RIKEN-BNL Research Center/BNL/ 
Columbia (RBC) Collaboration has recently fin- 
ished extensive quenched simulations of the 
hadron spectrum and low-energy QCD physics, 
using domain wall fermions || , with a particular 
focus on the residual chiral symmetry breaking 
effects due to finite values for the fifth dimen- 
sion, L s Q. (Similar work was also reported in 
0.) This work produced quantitative values for 
the residual quark mass m rcs , which enters low 
energy QCD physics through the effective quark 
mass, m c ff = wif + m Tes , where m/ is the in- 
put bare quark mass. The residual quark mass 
was shown to be small, ~ 1 /30 the strange quark 
mass for (3 = 6.0 and L s = 16. 

It was also found that topological near-zero 
modes of the domain wall fermion operator, which 
are not suppressed by a fermionic determinant in 
quenched simulations, have pronounced effects on 
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quark propagators for small to/. In addition, for 
large volumes where topological near-zero mode 
effects are suppressed, a non-linear dependence of 
to 2 on to/ was found, consistent with a quenched 
chiral logarithm. Thus, even though the chiral 
properties of domain wall fermions are quite im- 
proved over other formulations, these quenched 
pathologies can appear in simulation results. 

We have also finished measuring K — > it and 
K — > matrix elements for the 10 AS — 1 ef- 
fective four-fermion operators || relevant to de- 
terminations of e'/e and the real parts of the 
AI — 1/2 and 3/2 amplitudes for K ^ tttt de- 
cays. We have results for 200 quenched lattices 
of size 16 x 32 with L s = 16 and j3 = 6.0 and 
a variety of light and heavy quark masses. (For 
more details about the simulations see The 
analysis of these results requires control over the 
residual chiral symmetry breaking effects and an 
understanding of the role of quenched patholo- 
gies in our data. This report is a summary of our 
progress to date. 

A final step in determining physical matrix el- 
ements is the matching of continuum and lat- 
tice operators. This has been done using non- 
perturbative renormalization [[[o| and details can 
be found in fll|| . 

3. WARD IDENTITY CHECK 

Furman and Shamir || defined axial transfor- 
mations for domain wall fermions, which lead to 
the Ward-Takahashi identity 

A^(A«(x)0(y)) = 2m f (J«(x)0(y)) (1) 
+ 2(J« q (x)0(y))+i(8 a O(y)). 

Here A a ^(x) at the four-dimensional point x in- 
cludes spinors from all values of the the fifth co- 
ordinate. The pseudoscalar density J£(x) is con- 
structed from four-dimensional fields in the stan- 
dard way and J§ q (x) is constructed from "mid- 
point" spinors at s = L s /2 — 1 and L s /2. (The 
notation is as in Q.) 

For physics that involves momentum scales 
much less than the cutoff, the extra term in Eq. |l| 
involving Jg q (the "mid-point" term), gives a con- 
tribution of the form — m ICS Jg. This defines 
TO res to 0(a 2 ), since the low-energy QCD physics 



of domain wall fermions should by describable by 
an effective Lagrangian. However, in matrix ele- 
ments involving divergent quantities one only ex- 
pects that Jg q ~ 0(m les )J^. Since divergent ex- 
pressions enter in K — > -k amplitudes and we must 
manipulate the Ward-Takahashi identities to re- 
move them, these effects need to be understood. 

To test our understanding of TO rcs effects and 
our data, we first consider 0(y) = Jg(y) in Eq. 
|l|. Summing over x gives (no sum on a) 

]T([to/ JI{x) + Jg q (x)]Jg(y)) - (qq) = (2) 

x 

The term in the sum without J| ? gives to// to 2 
pion pole terms plus 0(m//a 2 ) contact terms 
along with constant terms and terms of higher 
order in to/. Similarly the mid-point term gives 
TOrcs/m 2 pion pole terms, plus 0(m YCS /a 2 ) con- 
tact terms as well as constant contributions. The 
pion pole terms enter as (to/ + m rcs )/m 2 in 
Eq. H while the contact contributions enter as 
0(to/) + 0(m res ). This form for the pion pole 
term is consistent with to w — at to/ = — TO ros . 

Thus we are lead to expect a Gell-Mann- 
Oakes-Renner relation for domain wall fermions 
with the form 

7)7 T n = - ??m, + 0-f +0 — 2-)(3 

2(TO/-|-TO ros ) a z a z 

This relation holds at non-zero values for to/, 
with the 0(m //a 2 ) contact term contribution be- 
ing cancelled by the similar term in {qq). Of 
course, there may be other quenched patholo- 
gies entering this expression, such as effects from 
topological near-zero modes or quenched chiral 
logarithms. 

We have measured m w , /, TO rcs and —(qq) for 
16 3 x 32 lattices at (3 = 6.0 J| and plot our re- 
sults in Figure [j]. One sees that for L s — 16, 
the larger mass points are approximately linear, 
but do not seem to extrapolate towards the value 
of —(qq) at to/ = — to ics . For L s = 24, the ex- 
trapolation is much closer, as is expected since 
TO r es has decreased by almost a factor of 2. Also 
shown are solid lines representing the left-hand 
side of Eq. |^ when a quenched chiral logarithm 
fit is used for to 2 (to,/). One sees that the ef- 
fect of chiral logarithms is to exacerbate the dis- 
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Figure 1. The GMOR relation for L s = 16 (up- 
per panel) and L s = 24 (lower panel). The open 
symbols are the left side of Eq. || and the filled 
symbols are — (qq). The dashed line gives the 
m / /a 2 dependence of — (qq) and the solid line in- 
cludes the effects of quenched chiral logarithms in 
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agreement in the extrapolation. Thus we have 
evidence for residual mass contact terms in this 
simple Ward-Takahashi identity. 

In conclusion, we note that for L s = 16, the 
quantity m 2 f 2 /(—2(qq)(mf + m rcs )), which is 1 
in chiral perturbation theory, differs from 1 by 
around 20% for our data. Thus one should avoid 
introducing extraneous powers of this quantity 
while manipulating the data. We now turn to the 
Ward-Takahashi identities useful for subtracting 
the divergent contributions to K — ► 7r matrix el- 
ements. 



4. TOWARDS 



\O n \K°) 



Here we discuss our tests of a particular method 
for removing the divergent effects in K — > 7r ma- 
trix elements. Following [fL2|, we define a sub- 



tracted operator 

O n = O n + rj n [{m d + m s )sd + (m d - m s )sj 5 d](A) 

where O n is a AS = 1 four-fermion operator and 
r]n is chosen such that 

(0\d n \K°) = 0. (5) 

Then following || , but noting that contact terms 
in the Ward-Takahashi identities coming from 
m ros effects will enter, we find 



(7r+7r-|0„|X°) = 



i(m 2 K0 -ml + )d{iT+\dn\K+) 



f 



dm 2 M 



(6) 



where the derivative removes any contact terms 
proportional to m Tes . 

The successful removal of the effects of mix- 
ing with the dimension three operators imple- 
mented in Eq. ^ requires chiral symmetry. We 
first check this with a simpler Ward-Takahashi 
identity, obtained by inserting an operator O = 
s(y)(l - J5)d(y) u(z)j 5 s(z) in Eq. |l|. The pion 
pole terms in the Ward-Takahashi identity then 
lead to 

(n\sd\K) 2 ^+r^ 

= l + 0(m f ) + 0(m ICS ) + --- 



(7) 



There are no power divergent operators in this 
Ward-Takahashi identity, so the coefficient of the 
rrif and m rcs contact terms are finite, but there 
may still be quenched chiral logarithm effects 
present. In Figure || we plot the left-hand side 
of Eq. ^ versus m/. We see that for larger m/, 
the points appear to extrapolate reasonably well 
to 1, but there is noticeable curvature for small 
quark masses. The figure includes results for 
axial vector sources, which agree with those for 
pseudoscalar sources, indicating that topological 
near-zero modes are not the origin of this non- 
linearity. We are investigating whether this can 
be explained as a quenched chiral logarithm. 

We determine r\ n from measurements of 
(0\O n \K°) for non-degenerate masses for the 
strange and down quarks. The upper panel of 
Figure | shows a plot of (0|O 2 |^ }/(0|s7 5 d|iY } 
versus m s — md- We see good linearity which al- 
lows for a reliable determination of 772 , with sta- 
tistically errors at the few percent level for this 
particular operator. 
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Figure 2. The left-hand side of Eq. |?] plotted 
versus quark mass. 



With the value of 772 determined, 2 is known 
and Figure |^ shows 
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rj 2 [m s + m d/ 



sd\K+) 



(tt\P\0)(0\P\K) 



(8) 



where OJ is the Al =1/2 part of 2 - In lowest 
order chiral perturbation theory, the denomina- 
tor is a constant, so the desired K — > tttt matrix 
element is given by the slope and the intercept 
is related to m rcs . We are looking carefully for 
non-linearities in our data, paying particular at- 
tention to the quenched chiral logarithm effects 
recently calculated in Jig] ]. 

5. CONCLUSIONS 

We have seen that the Ward-Takahashi identi- 
ties for domain wall fermions include 0(m re s/a 2 ) 
contact terms which can be removed by taking 
derivatives with respect to mj. Our data also 
show evidence for non-linear dependence on mj, 
some of which may be due to quenched chiral log- 
arithms. We have shown that the subtracted op- 
erator 2 can be determined and are now doing 
the subtractions for all the operators. 

The calculations reported here were run on the 
QCDSP computers at Columbia University and 
the RIKEN-BNL Research Center. 
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Figure 3. The value for (0\O 2 \K°) / {0\sj 5 d\K } 
is shown in the upper panel. The lower panel 
is a plot of the K — > 
normalized as in Eq. ^ 
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